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A Submillimetre Imaging Polarimeter at the James Clerk Maxwell 
Telescope 
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ABSTRACT 

A polarimeter has been built for use with the Submillimetre Common-User Bolometer Array 
(SCUBA), on the James Clerk Maxwell Telescope (JCMT) in Hawaii. SCUBA is the first of 
a new generation of highly sensitive submillimetre cameras, and the UK/Japan Polarimeter 
adds a polarimetric imaging/photometry capability in the wavelength range 350 to 2000 j^m. 
Early science results range from measuring the synchrotron polarization of the black hole 
candidate Sgr A* to mapping magnetic fields inferred from polarized dust emission in Galactic 
star-forming clouds. We describe the instrument design, performance, observing techniques 
and data reduction processes, along with an assessment of the current and future scientific 
capability. 

Key words: 

instrumentation: polarimeters — techniques: polarimetric — submillimetre — polarization — 
magnetic fields 



1 INTRODUCTION 

Linear polarization of submillimetre continuum radiation is as- 
sociated with optically-thin synchrotron emission, and with par- 
tial magnetic alignment of elongated dust grains. In both cases, 
magnetic morphology can be deduced from the directions of po- 
larization vectors, and these techniques are much less subject to 
confusion than in the optical and radio regimes, where scattering 
and Faraday rotation dominate the polarized signal. In the submil- 
limetre, the limitations are that the data are only sensitive to the 
net plane-of-the-sky magnetic field direction within the telescope 
beam, and that there is no direct information on magnetic field 
strengths. However, with the capability to detect levels of linear 



polarization of 1% or less, considerable magnetic structure infor- 
mation can be obtained. The importance of magnetic fields is now 
being realised in sources ranging from disks around young stel- 
lar objects to the inner jet s in active galactic nuclei iTamura et alJ 
ll999tlNartallo etalll998l) . 

Observing techniques at submillimetre wavelengths have ma- 
tured considerably in the last decade, particularly with the intro- 
duction of the first generation of continuum cameras. The chief dif- 
ficulty for ground-based observations is the high opacity and time 
variability of the Earth's atmosphere, and these problems are crit- 
ical for polarimetry. Atmospheric absorption can be largely over- 
come by building telescopes at suitable high, dry sites, such as 



2 J. S. Greaves et al. 



e~ Chopper Drive Motor 
— Chopper Wheel 

r— Chopper Control Box 



Rotating WavQplate 
ever SCUBA Window 




Drive Pulley It Motor 



Polarimeter Casing 
Encoder Wheel 



Chopper Wheei 

Chopper Control Box 




Insulated Mounting Plate 
Wedged Polarteer/Analyser 
Waveplate Assembly 

16 Position Encoder Disc 
Opto Sensors 



-Cryostot 
Scuba Window 
Scuba Filter Drum 



Figure 1. Technical drawing of polarimeter (from Murray et all ll997t) '). The physical dimensions of the polarimeter module are 30 by 10 cm (maximum 
length and thickness respectively). The waveplates are 12 cm in diameter; the thickness of individual half-wave quartz plates depends on wavelength but a 
representative value is 0.9 cm for a single plate optimised for 850 ytm observations and made of quartz with a submillimetre refractive index of 2.1. 



Mauna Kea in Hawaii (summit altitude 4200m), and by observing 
in frequency 'windows' away from oxygen and water absorption 
lines. Under the best conditions, the zenith transmission is then 
about 80% at a wavelength of 850 iim, and as much as 40% at 
shorter wavelengths such as 450 \xm. The other problem affect- 
ing the quality of flux measurements is atmospheric instability, es- 
pecially emissivity changes on short timescales of a few seconds. 
This is a particular challenge for polarimetric observations, which 
require measurements accurate to better than 1% over periods of 
minutes to hours, and necessitates sophisticated techniques for re- 
moving sky-level changes. 

Submillimetre polarimetry using single-pixel photometers be- 
gan about a decade ago (e.g.. lFlett & Murra^ll99lh . and a num- 
ber of important detec t ions were made both in th e (sub)millimetre 
(e.g jKane et all 1 9931: iGreaves et alll997llGlenn et alll999h and 
with complementary studies of warmer sources in the far-infrared 
(e.g.. lNovak et alJl989» . These observations were extremely labo- 
rious, but a major breakthrough has been made with the arrival of 
sensitive submillimetre and far-infrared cameras. With the inclu- 
sion of polarimeters — typically based on rotating half-wave plate 
designs and either internal or external to the cryogenically cooled 
instrument — polarization imaging of a wide variety of sources 
has now become feasible. A complete introduction has been pre- 
sented by Hilde brand et alj fcOOd) . In this paper, we describe ob- 
servations and techniques with the new imaging polarimeter 1 used 

1 The SCUBA Polarimeter has been funded jointly by the UK Particle 



with SCUBA, the Subm illimetre Common-User Bolometer Ar- 
ray jHolland etaD ll999) at the James Clerk Maxwell Telescope 
(JCMT) on Mauna Kea. 



2 INSTRUMENT DESIGN AND PERFORMANCE 

The general design of the polarimeter has been described by 
iMurrav et alj ll997l) . It is an external module that fits over the 
SCUBA cryostat window, and comprises a rotating half-waveplate 
and a fixed photo-lithographic grid ('analyser') (Figure 0. The 
waveplate is stepped to a series of fixed angles, which has the effect 
of rotating the incoming source plane of polarization, by twice the 
physical angle of rotation. The grid, with a 10 micron line spac- 
ing (of etched copper on a mylar substrate), efficiently reflects one 
plane of polarization that is fixed in the SCUBA reference frame. 
The net effect (Figure|2j is to select out a fixed component of a ro- 
tated linearly polarized signal, so that the SCUBA bolometers see 
a modulated signal of the form 

S(S) = l/2[/„ + I p (l + cos(4S - 26))]. (1) 



Physics and Astronomy Research Council and the NAOJ and JSPS of Japan. 
SCUBA was funded by the JCMT Development Fund supported by PPARC, 
the National Research Council of Canada and the Netherlands Organisation 
for Pure Research, and built by the Royal Observatory Edinburgh, now the 
UK Astronomy Technology Centre. 
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Figure 2. Illustration of the effects of a rotating half-waveplate and analyser 
iBerrv & Gledhill200fj|) . 



Table 1. Waveplate properties (symbols defined in the text). The input sig- 
nal was produced with a second analyser to make a pure linear polarization 
from an astronomical source (generally Saturn). The PME is a function of 
waveplate angle with respect to the incoming plane, but minimum and mean 
values differ by no more than 8%. The <j> values in brackets refer to obser- 
vations after late 1998, when the photometric plate was re-mounted; the 
1350 (im value was measured and the shift at 2000 |am was inferred. 
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Here I u is the unpolarized intensity, I p describes the modulating 
polarized component above the minumum observed signal I u , 8 is 
the waveplate angle, and 6 is the position angle of the polarization 
in a suitable co-ordinate frame. The factor of 1/2 arises because 
one plane of polarization is reflected by the analyser; this contrasts 
with cameras designed specifically for polarimetry, which typically 
use a waveplate inside the instrument and measure both planes 
of po larization using two separate detector arrays jPowell et alj 

2.1 Waveplate design 

The waveplates are made of birefringent quartz, cut to a thick- 
ness to retard one plane of polarization by half a wavelength rel- 
ative to the orthogonal plane, on emerging from the plate. Since 
SCUBA has filters for observing at a wide range of wavelengths, 
from 0.35 to 2 mm, an ac hromatic design was adopted foll owing 
the Pancharatnam method iMurrav et al Jll997tlritie" & Rosenberg 
Il98ll) . This uses a sandwich of an odd number of individual half- 
waveplates whose 'fast' axes of polarization are offset by ±60°. 
The effect is a combined plate whose fast-axis direction (</>) is 
slightly wavelength-dependent, but with very good polarization 
modulation efficiency (PME) over a broad range of wavelengths. 
The plates were wax-bonded and then held in place in mounting 
rings. A good transmission (77) was achieved using single-layer 
polypropylene coatings to reduce reflections at the plate surfaces. 

For the SCUBA polarimeter, two achromatic designs were 
used: one with three quartz plates to cover the 1.1, 1.35 and 2 mm 
bands and one with five plates for the 350, 450, 750 and 850 \im 
bands. Because SCUBA has single photometers at the millimetre 
wavelengths and full imaging arrays at the submillimetre wave- 
lengths, the plates are known as 'photometric' and 'array' for con- 
venience. 



not be removed to change the plate.) Nominally, measured polariza- 
tion percentages should be divided by the PME to obtain the true 
values, but in practice this small correction is usually neglected. 
The transmission properties met the specifications for a single-layer 
coating, except at the shortest observing wavelengths, and it is sus- 
pected that re-coating and bonding of the plate might be necessary 
to improve this. 

2.3 Reproducibility 

An important part of the commissioning test was to confirm that 
the results agreed with previous polarimetric observations, from the 
JCMT and elsewhere. Tests included both low and high polariza- 
tion sources: for example, 

• 850 |i.m polarization of the DR21 cloud core was measured to 
be 1.70 ± 0.27% at 22 ± 4°, whereas Minchin & Murray (1994) 
found 1.8 ± 0.3% at 17 ± 4° at 800 \im 

• a flux peak in the Crab Nebula was de tected at 1350 \im with 
25.9 ± 0.6% polarization at 148 ± 1 . while lFiett & Murray! \l 99 ll) 
obtained 25 ± 2% at 144 ± 2° at 1 100 urn. 

These results were from previous photometric polarimetry at the 
JCMT, and are for dust and synchrotron emission respectively. 
Since the previous instrument had no imaging capability, repro- 
ducibility in this mode relied on data from other telescopes: for 
exam ple, our 850 \im polarization map of OMC1 1 CooDin et al. 
2000 ) is very similar to the 350 |J.m results of ISchleunind 1 1998, 
his Figure 5). The position angles of dust polarization should be 
wavelength-independent provided that the telescope beams sam- 
ple th e same grain population and line-of-sight 1 Hildebra nd et alJ 



2.2 Performance on the JCMT 

The performance of the waveplates on the telescope was first mea- 
sured in October 1997, and the results are given in Table 1. Excel- 
lent PME values were measured, close to 100% at all wavelengths, 
and it was even found possible to operate the photometry plate in 
second order (3 A/2) at the array wavelengths. (This has advantages 
for complex observing programmes, because the polarimeter need 



2.4 Sensitivity limits 

The sensitivity for polarization measurements is determined essen- 
tially by the sensitivity of SCUBA, with additional effects of the 
transmission through the polarimeter and the limiting accuracy of 
corrections for sky fluctuations. The critical quantity is the polar- 
ized flux: the product of percentage polarization and mean flux. 
With the previous JCMT polarimeter used on a single-bolometer 
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Table 2. Percentage and position angle (p, 8) of the main-beam instrumental polarization (not corrected for PME or tj>); for the array data the uncertainties 
are the standard error of the set of individual measurements. Observations were made using several planets; note that Saturn may be slightly polarized due to 
scattering by the rings (with an best-fit value of 0.5% at 850 urn); however this effect is greatly reduced in data averaged over many hour angles. Saturn is also 
larger than the beam size at all the submillimetre wavelengths, so these data include some degree of sidelobe polarization. 



A 



p(IP) 

(%) 



0o(IP) 
(degrees) 



epoch 



Array waveplate ( central array bolometer) 

350 0.82 ±0.13 100 ±4 

450 3.51 ±0.34 113 ±3 

750 1.16 ±0.04 97 ±1 

850 1.09 ±0.06 161 ±2 
Array waveplate (off-centre bolometers) 

450 3.26 ± 0.23 100 ± 2 

850 0.92 ± 0.05 163 ± 2 
Photometric waveplate 

1350 1.72 ±0.11 166 ±2 Uranus 
2000 1.34 ±0.09 170 ±1 Saturn/Uranus 



Saturn 
Mars/Uranus 

Saturn 
Mars/Uranus 



Saturn 
Mars/Uranus/Saturn 



Oct 1997 + Aug 1999 
1997-1999 (5 dates) 
Oct 1997 + Aug 1999 
1997-1999 (25 dates) 

Jul 1998 (6 bolometers only) 
1997-1999 (4 dates, all 36 bolometers) 

1997-1999 (7 dates) 
Oct 1997 & Aug 1999 



detector, UKT14, the smallest polarized flux that could be mea- 
sured was about 100 mjy at 800 \im iMurrav et al Jl 19971) . This is 
equivalent, for example, to p = 1% in a 10 Jy source, such as one of 
the brightest Galactic protostars. 

With the SCUBA polarimeter, the expected sensitivity can be 
derived by considering the magnitude of the fractional polariza- 
tion, as obtained by subtracting signals at two waveplate angles 
and dividing by their sum 2 . The error on the polarization fraction 
is 1/ v2 times the fractional intensity error, so half as much time 
is needed as would be for a measurement of I alone. However, this 
is cancelled by a factor of twice as much time that arises because 
half the photons are lost to reflection off the analyser. Finally there 
is another time factor of four because four waveplate positions are 
needed to obtain the two orthogonal linear polarization components 
(Stokes parameters). This determines the signal-to-noise on the po- 
larization percentage, <t p , for a given integration time t: 

t = 4/r/ x N x (NEFDtjp/pF) 2 (2) 

where NEFD and F are the noise equivalent flux density and source 
flux in mJy/\/Hz and mjy respectively, and p is the fractional polar- 
ization (e.g. 0.01 for 1%). The factor N is the inverse of the fraction 
of time spent on one spatial point, so it is 1 for photometric po- 
larimetry, 4 for imaging polarimetry at beam width spacings, and 
16 for a Nyquist-sampled polarization map. 3 

At the primary observing wavelength of 850 p.m, the NEFD 
is 70 mJy/\/Hz under the best conditions, and a 1 hour integration 
would give a 3a detection for a polarized flux p x F = 15 mJy in a 



2 Two waveplate angles 45° apart will have the largest and smallest ob- 
served signals, I m ax and I m in, and it can be demonstrated from Equa- 
tionQthat {Imax ~ Imin) /(Imax + Imin) is equivalent to the polar- 
ization fraction, p = I P /(I U ± Ip)- For convenience we re- write p = 
(R. — 1)/(R + 1) where R = Imax / Imin, and differentiate this to give 
dp = 2/(R + l) 2 dR. For small polarizations R 1 so dp dR/2. Then 
as the two / measurements are independent, dR/R = \/2dI /I and com- 
bining the last two expressions, dp = dl/I X 1/ \/2 in the limit H — * 1. 

3 The SCUBA beams are spaced two beam widths apart on the sky, so that 
at beam width resolution a bolometer is effectively looking at a particular 
point for 1/4 of the time. Since the array is undersampled, the telescope 
secondary mirror is 'jiggled' to fill in the image; the polarization maps are 
usually sampled every 6" but smoothed to about beam width spacing to im- 
prove o-p. For an extended source, smoothing with e.g. a Gaussian function 
includes more flux and hence reduces t to obtain a given a p . 



beam-spaced map, or 7.5 mJy for photometry. This signal-to-noise 
is sufficient to determine magnetic field directions to 0.5 radians 
/ 3 jNashizadeh-Khouei & Clarkdll993l) . or about ±10°, and the 
sensitivity is more than an order of magnitude better than with the 
UKT14 polarimeter. The faintest sources actually detected have had 
fluxes of about 0.2 Jy per beam and polarizations of a few percent; 
scientifically, large samples of protostars, faint Galactic clouds and 
AGN can now be detected. Also, types of sources that were pre- 
viously impossible can now be studied, including asteroids, pre- 
stellar cores and nearby starburst galaxies. 



3 INSTRUMENTAL LIMITATIONS 

3.1 Main-beam instrumental polarization 

The absolute accuracy of any polarization measurement depends 
on how reliably instrumental effects can be subtracted. For the 
SCUBA polarimeter the main problem is instrumental polarization 
(IP), which is dominated by the woven Goretex windblind through 
which the JCMT observes. The thread spacing, which is approxi- 
mately 0.5 mm, is slightly different in the vertical and horizontal di- 
rections and this affects the relative transmissions, particularly at a 
wavelength of 450 \im which is closest in size. The solution would 
be to either roll back the windblind (a half-hour labour-intensive 
operation that is only possible on rare occasions of very low wind- 
speed) or add a second orthogonal sheet of Goretex in the beam to 
induce IP cancellation. The latter has been tried and reduced the 
IP by a factor of at least three at 450 \im; further experiments are 
ongoing but this is not yet the standard observing mode. 

The IP in the main beam of the telescope has been mea- 
sured using planets, assumed to be unpolarized (e.g. Jclemens et all 
1990), and the results are listed in Table 2. At the primary wave- 
length of 850 \im the measured IP values are very stable and can 
be subtracted to an accuracy of about ±0.25%, the measurement 
uncertainty for individual bolometers. Global differences, such as 
off-axis effects in the three rings of bolometers that form the hexag- 
onal close-packed array, are 0.1% or less, and for the central array 
bolometer (used for photometric polarimetry) the measurement er- 
ror has been reduced to about ±0.06%. Fewer data are available at 
other wavelengths but generally a very accurate measurement can 
be made during a particular observing run. Values tend to be sta- 
ble over time, taking into account that the angles must be corrected 
for telescope elevation. The observed 9 is So plus elevation, due to 
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Figure 3. Typical 850 |am map of polarization measured around Jupiter, 
due to sidelobe IP. Vectors shown range from 0.8% to 15% (with criteria 
p/Sp > 2 and 8p < 3%) and are plotted every 6" (2 pixels). Main-beam 
instrumental polarization has been subtracted. Data are from January 2000, 
using a 1 80 arcsec chop of the secondary mirror in the azimuth direction. 

the alt-az telescope mounting and the fact that SCUBA is aligned 
with the elevation axis but is fixed on the Nasmyth platform. This 
angle relation is quite reliable, with deviations of around ±15° at 
extreme elevations, for the 850 \im pixels at the array edges that 
are furthest from the optical axis. 



3.2 Sidelobe polarization 

The sidelobes of the telescope beam have low intensity-amplitudes 
(e.g. < 1% at 850 [im at > 40" from Jupiter in Figure|5J, but differ- 
ent reflection and scattering effects for the two incoming planes of 
polarization can still produce instrumental 'sidelobe polarization'. 
This needs to be considered when observing extended sources, or 
faint objects next to much brighter ones, as the sidelobe IP can be 
considerable, and is additional to the main-lobe IP that has been 
measured using a compact planet. 

In theory, a polarization map of an extended unpolarized 
source (generally Jupiter), made with the same primary dish sur- 
face and secondary mirror chop throw as the source observations, 
can be used for an exact subtraction of sidelobe IP. An example 
of such a map is shown in Figure [3] after subtraction of the stan- 
dard main-lobe instrumental polarization. However, in practice it is 
often easier simply to estimate the error induced by the sidelobe ef- 
fects. The most extreme case is seen with large chop throws of the 
secondary mirror (the maximum reccomended is ±180") in a non- 
azimuthal direction. In this case, there are residual vectors towards 
about half of Jupiter's disk, with magnitudes up to 3%, although 
when averaging over this area (diameter of 40 arcsec) the net effect 
is only 0.25%. For smaller or azimuthal chop throws, the net resid- 
ual is similar or lower and only about 1 in 10 individual points has 
a significant vector. Thus if the observer is interested in average po- 
larization of an extended source, the error is quite small, but more 
careful calculations are required for data on a single map point in 
an extended region. 



Additional problems arise if trying to extract source polar- 
ization when there is a much brighter object nearby in the field. 
In Figure|5] there are vectors averaging ~ 6% around Jupiter, al- 
though the largest tend to be at large radii where the power in the 
beam is low. A good general method is to set a threshold of source 
polarization percentage, above which sidelobe IP is not a serious 
contaminant, and to base this threshold on the relative polarized 
flux contributions. For example, if the instrumental polarized flux 
is half that of the source polarized flux, then in the worst case the 
two effects have orthogonal Stokes parameters, and the error in 9 
for the source will be 1/2 tan -1 0.5 or only 13°. (If the effects are 
parallel, then the error will instead be in p, which is generally of 
less scientific importance.) 

Taking into account the beam power profile, this criterion 
for the minimum believable source polarization percentage can be 
written 

p crl t>2x Ps ,P s i(F s i/F) (3) 

where the three terms on the right-hand side are the IP at the rele- 
vant point in the sidelobe, the power here relative to the main beam, 
and the ratio of the flux of the object in the sidelobe to the flux in 
the main lobe. Thus sidelobe polarization can often be non-critical 
if the first two terms are small, even if there is a bright object 
in the sidelobe (F 3 i >> F). As a quantitative example. iGreavesI 
(2002) observed polarization in the jets of the proto-planetary neb- 
ula CR 2688, at positions offset 15" from the central star. The ratio 
F s i /F was 3.6, the average p s i was 3% and the power at the edge of 
the main beam was P s i w 0.06. Using Equation|3| the polarization 
of the southern jet was found to be believable because the measured 
value of 1.9% exceeded p cr it evaluated at 1.3%. At larger offsets, 
the critical product of power and IP tends to be smaller: for exam- 
ple the values < 0.01 and ^ 15% (Figure[3j give a smaller product 
than in the CRL 2688 case. The most critical quantity in individ- 
ual sources is thus the flux ratio: very high contrasts will mean that 
faint regions need to be highly polarized to give usable results. 

4 OBSERVING TECHNIQUES AND DATA REDUCTION 

Data acquisition and reduction are based on extensions to standard 
techniques used at the JCM T; more details of the stan dard proce- 
dures are given in Holl and et all <1998Ul999l) . ljenness & Lightfootl 
( 1998) and Jenn ess^^lTl20001) . For polarimetry, standard imag- 
ing or photometry observations are made at a set of waveplate po- 
sitions (normally 0, 22.5, 45, 67.5, 90... 337.5°), and these data are 
fitted to detect the sinusoidal signal modulation from rotating the 
incoming polarization vector. Using a complete cycle is intended 
to eliminate other modulations such as reflections off the wave- 
plate. An alternative to fitting is direct subtraction of data taken 
at different angles, as for example, p cos(26) can be found from 
(5(45) - 5(0))/(5(45) +5(0)) andp sin{26) from (5(67.5) - 
5(22.5))/(5(67.5) + 5(22.5)) (Equation [TJ; this method is en- 
tirely valid and gives both Stokes parameters directly, but is less 
often used. For an excellent introdu ction to polarization observing 
procedures and limitations, see also Hildebrand et al. (2000). 

4.1 Observing modes 

Three observing modes have been adapted for use with the po- 
larimeter. These are photometry, where a single bolometer observes 
a fixed point; jiggle-mapping, where the secondary mirror is jig- 
gled to obtain an image of one SCUBA field-of-view, using all 
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the bolometers in an array; and scan-mapping, where the telescope 
scans across a source and the data from all bolometers are recon- 
structed to give large images. All of these modes use chopping of 
the secondary mirror to remove the sky emission; photometry and 
jiggle-mapping also nod the source between left and right beams 
halfway through an observation to minimise sky gradients. For po- 
larimetry, the only difference is that a complete standard observa- 
tion is made at each waveplate angle before moving on to the next 
angle. 

The first two polarimetry modes are fully commissioned, and 
the third has been proven experimentally but is not yet fully in- 
tegrated. All six of the available SCUBA wavelengths have been 
used, but only a small number of observations have been made with 
the 350/750 ixm filter combination or using the 2 mm photometric 
pixel (less than 1% each of the total number of observations). Thus 
the best quantified performance is at 850 and 1350 450 |i.m 
data are obtained simultaneously with 850 \im but are often not us- 
able due to poor atmospheric transmission. Scan-mapping has been 
tested only at 850 \im. 

The photometry mode is most suited to sources with accu- 
rately known positions that are smaller than the telescope beam 
width (8, 15 and 24 arcsec at 450, 850 and 1350 \im). The off- 
source bolometers remain 'switched-on' and can be used to subtract 
residual sky signals (not possible with the individual photometric 
pixels at 1350 and 2000 |J.m). Typically the secondary mirror chops 
between on- and off-source with a 1 arcmin throw, to allow sub- 
traction of the sky DC-level, and 8 seconds of data are obtained at 
each of the 16 waveplate positions. The chop takes place at a rate 
of 7.8 Hz, and the telescope nods so that the source moves from the 
left beam to the right beam after 4 seconds. A small grid of points 
can be used rather than a fixed point but this in not generally done. 
Integrations of up to about 20 seconds can provide higher signal- 
to-noise in stable conditions, but generally many short integrations 
are co-added. 

The jiggle map mode is typically used for sources bigger than 
the beam but smaller than the 2.3 arcmin array field-of-view in at 
least one dimension. If the source is larg er, two effects can cor- 
rupt the polarization measurement (Matthews et al. 2001). Firstly, 
the maximum recommended chop throw is 3 arcmin, so sources 
larger than this will be 'self-chopped'. In practice sky noise sub- 
traction requires several emission-free bolometers, so the setting of 
the zero level will be inaccurate if the source dimensions approach 
3 arcmin; an incorrect zero-level changes I u and so alters the per- 
centage polarization. Secondly, if there is emission in the reference, 
i.e. nominally off-source, beams it may well be polarized and this 
can corrupt the on-source polarized signal in both magnitude and 
direction. A limit can be estimated where the reference polarization 
is not known, but this may not set useful constraints: for example, 
faint outer regio ns in some dust clouds can be up to ~ 15% polar- 
ized (e.g jMatthews et all200ll) . 

A 16-point jiggle sequence with the secondary mirror gives an 
850 (or 750) vim image with a point every 6 arcsec, so is slightly 
better than Nyquist (half-beam) sampled. The simultaneous data 
at 450 (or 350) \im are undersampled, since the points are only a 
little better than a beam width apart, and a 64-point jiggle sequence 
with 3 arcsec sampling is required if these data are of importance. 
The simpler jiggle takes 32 seconds to complete at each waveplate 
position, and the more complete jiggle takes 128 seconds, so the 
former is preferable for greater sky stability. 

The scan map mode has been tested at 850 \im only, on a few 
sources with sizes of about 6-10 arcmin. The technique used is to 
scan the telescope across the source in any suitable direction, while 



chopping in RA or Dec. Using two different chop throws in each 
direction, for a total set of four, has been found to give good recon- 
structed images; further det ails of this 'Emerson IF technique are 
given by Jenness et al. 1 2000). It is important to scan off-source so 
that a true baseline can be established and set to zero in each indi- 
vidual waveplate image, before the polarization reduction is done. 

The scan rate used is 48 arcsec per second, so sampling ev- 
ery l/8th second gives points 6 arcsec apart, the same as in jig- 
gle mapping. The minimum number of waveplate positions needed 
to deduce the Stokes parameters is four, so scan map polarimetry 
is rather time-consuming, for example, taking 30 minutes for one 
dataset on a 400 arcsec square area. This compares to three minutes 
or less for a set of Stokes parameters in the other modes, so scan 
map polarization measurements require more stable conditions. 



4.2 Data reduction techniques 

Two separate software packages are used for the da ta reduction: 
one for photometric polarimetry (SIT:lNar talldll9 95h and one for 
imaging polarimetry (POLPACK: lBerrv & Gledhilll2000h . The re- 
duction philosophy is the same in both cases, and in fact POLPACK 
can handle single-pixel 'images', while SIT is more optimised for 
inspecting the results of individual waveplate cycles on a single 
point. 

The raw data are reduced Ijenness et alj|2002l) by extracting 
the chopped signals, flatfielding, and correcting for atmospheric 
extinction using a skydip. This is a series of measurements of 
sky emission at different elevations, made either with the JCMT 
at the observing wavelength, or taken from the 1.3mm database 
recorded at the adjacent Caltech Submillimeter Observatory (CSO) 
and extrapolated in wavelength. The next refinement is to subtract 
rapid changes in sky emissi on (to a greater accuracy than achieved 
simply by chopping: l Holl and et alj n"999)) by shifting the mean 
level using blank bolometers in the array, or in scan-mapping, sub- 
tracting smoothed data seen by each bolometer over 2 seconds. 
This 'sky noise' is spatially correl ated over the array diameter 
IJenness. Lightfoot & Holland 19981) . 

The main lobe instrumental polarization is then removed on a 
bolometer-by-bolometer basis using the expression: 



s 



measured 



l+pip{e) cos(4<5-26' / p(e)) 



(4) 



where pip(e),8ip(e) are the percentage and direction of instru- 
mental polarisation at elevation e (a constant percentage is gener- 
ally assumed). This is an approximation of 

Scorrected = Smsasured ~ S mea nPlp{e) COs(4<5 — 28[p(e)) (5) 

and is valid for instrumental polarizations of up to a few percent. 
The former expression is prefered as the mean flux level is not eas- 
ily determined prior to regridding — the bolometer jiggles to differ- 
ent source positions and thus only part of the data stream refers to a 
particular spatial point. Also, not exactly the same point is seen by 
this bolometer for the next waveplate angle, because of sky rotation 
(SCUBA does not use an image de-rotator). 

In imaging mode the data are then regridded to a rectangular 
array, removing pixels from the edge that are affected by regridding 
edge effects in order to simplify mosaicing. For scan map data, a 
Fourier deconvolution is used to remove the chop sig nature and 
combine the four different chopped maps Ijenness et all2000l) . 

The Stokes parameters I, Q, U are then extracted by fitting 
the signal in each pixel as a function of waveplate angle; errors 



are supplied by comparing signals at identical or equivalent wave- 
plate positions (e.g. 0, 90, 180 and 270°; see Eq.[TJ) In SIT, the 
Q and U are then listed and statistical tests can be used to reject 
anomalous fits before calculating averages and deducing p, 6. In 
POLPACK, an I, Q, U data-cube is generated and a catalogue of 
vectors produced corresponding to all the image pixels. When com- 
bining different data sets (a set contains 16 waveplate positions) it 
was determined that the best results were obtained by coadding the 
IQU cubes rather than adding more data to the fit. As more data 
were added to the fit the signal-to-noise did not increase in the ex- 
pected way, suggesting that it was susceptible to DC-level differ- 
ences amongst the observations. 

Various binning and selection procedures can be used to im- 
prove the polarization signal-to-noise. Typical criteria are p/5p > 
3 (i.e. 86 < 10°) and an upper limit such as 5p < 1 — 2% to 
eliminate biased data. Because the waveplate angles increase step 
by step, a linear drift in the sky transmission can force, for ex- 
ample, a positive Q and U, and these spurious fits tend to pro- 
duce a high p, especially in regions of low I. Thus a Sp criterion 
can eliminate such data with a large scatter between good and in- 
accurate fits. For example, four measurements with similar 9 but 
p — 3, 3, 3, 9% would give an average p of 4.5% and a standard 
error of 1.3% (hence p/Sp = 3.5). This dataset is clearly biased by 
the last measurement, but could be eliminated with a 5p < 1% 
criterion. This entire reduction process h as been automated us- 
ing the ORAC-DR data red uction pipeline lEconomou et all 19981 ; 
Ijenness & Economoull999T) . 



5 SCIENTIFIC CAPABILITY 

Science targets so far observed include asteroids, planetary 
nebulae, supernova remnants, T Tauri star disks, accreting 
protostars, pre-stellar cores, Bok globules, dark cloud fila- 
ments, high-mass star-forming clouds, the Galactic Centre in- 
cluding the black hole Sgr A*, the starburst galaxy M82, 
the ultraluminous galaxy Arp 220 and a number of vari- 
able AGN o r 'bla z ars'. Results have been di scussed by 
iTamura et all ll999h: iMatthews & WilsorJ 120021): iDavis et alj 
fcOOCh: IWard-Thompson et alJ feOOOh; [Henning et alJ | 200lT) ; 

[ Vallee et alJ <200d). iFeldman et alJ feOOOfl : lAitken et alJ kood) ; 
Greaves et al. 1 2000); Marscher et al. 1 1999), among others. Many 
of these objects have sub-Jy fluxes per beam, and were not feasi- 
ble to observe before SCUBA; they are also very difficult for other 
ground-based instruments such as the HERTZ polarimeter-camera 
used at the CSO, which operates in the more challenging 350 u.m 
band iDowell et all 199 8). Examples of some typical observations 
are described below; the imaging results shown in the figures have 
not previously been published. 

5.1 Polarimetric Photometry 

Sgr A*, the massive black hole candi date at the centre of the 
Galaxy, was observed polarimetrically by Aitken et al. 1 2000). Cru- 
cial observations were made at 1350 and 2000 u.m in the photomet- 
ric mode, to establish the polarization spectrum in the millimetre 
regime as well as at submillimetre wavelengths. These were the 

4 Scan map observations are sufficiently time-consuming that only the an- 
gles 5 = 0, 22.5, 45, 67.5° may be used, in which case the same angle in 
different cycles should be compared. 



Submillimetre Imaging Polarimeter 7 




Figure 4. Fits to 2000 pm photometric polarimetry of Sgr A*, illustrating 
data quality for one waveplate cycle on an 8 Jy, 3% polarized source (in- 
strumental polarization of 1.4% has not yet been subtracted). The axes are 
measured signal modulation versus waveplate orientation, and the scales are 
polarization fraction (tick marks at intervals of 0.01, i.e. 1%) versus angle 
(tick marks at 10° intervals). The top plot covers 5 = 0° — 157.5° and the 
bottom plot is for 8 = 180° — 337.5°. Good consistency is demonstrated 
by the upper and lower fits from the first and second halves of the waveplate 
cycle being in phase. 

first detections of linear polarization of Sgr A*, and they showed a 
dramatic ~ 90° shift in position angle around 1 mm wavelength, 
which implies self-absorption in a very compact emission region of 
only a few Schwarzschil d radii 1 Aitken et alJl200ollBromlev et alJ 
l200ltlLiu & Melij2002h . The 1350 and 2000 u.m detections were 
each significant at the 10 sigma level after about 15 minutes of inte- 
gration. An example of the fitted results from one waveplate cycle 
(2 minutes of data) is shown in Figure|4| 

It is now routinely possible to detect polarized emission from 
compact radio sources. A ~ 5a detection of a 0.5-10 Jy source can 
be made in 30-60 minutes, depending on the fractional polarization 
(values of 2-40% have been observed). Large levels of polarization 
indicate that the magnetic field must be quite highly ordered, and 
allow its direction to be inferred without worrying about opacity 
effects. Comparisons with quasi-simultaneous results from 7 mm 
very long base line interferometry show that, in many cases, the 
angle of submillimetre polarization matches that very close to the 
base of the sub parsec-scale jet, and/or the structural position angle 
of the jet (Stevens et al. in prep.). 

One of the faintest photometric detections with the 
UKT14 polarimeter was the HH24MMS protostar observed by 
Greaves, Holland & Ward-ThomDson ( 1997) at 800 [j.m. This re- 
sult was confirmed with the SCUBA polarimeter at 850 [im: 
HH24MMS was 3.7 ± 1.3% polarized at 95 ± 10° in the earlier 
data and the new measurement was 3.82 ± 0.61% at 109 ± 5°. 
Thus results published with the old instrument appear to be quite 
reliable, but now much fainter targets can be detected. 

5.2 Polarimetric Imaging of Single Fields 

Jiggle map polarimetry has been used for a wide variety of tar- 
gets: Figure [5] illustrates results for the high-mass star-formation 
region W3. The two bright cores contain the luminous clusters of 
young stellar objects IRS5 and IRS4 (left and right respectively), 
plus an extended envelope. The image illustrates the care that must 
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Figure 5. 850 |am polarization map of the main cores in W3 (see text). The 
colour scale goes up to 8 Jy/beam and the vectors range from 0.4% to 7.6%; 
co-ordinates are B1950. All detections have p/Sp > 3 and Sp < 1%. The 
i?-plane of polarization is plotted, i.e. parallel to the long axes of aligned 
dust grains and perpendicular to the aligning magnetic field. 




R.A, offset (arcsec) 

Figure 6. 850 |am polarization data for the Crab Nebula, using scan map 
polarimetry. The smallest and largest polarization vectors are 8% and 58%, 
and the colour scale goes up to 1.8 Jy/beam at 850 \xm. There are 183 vec- 
tors with p/Sp > 5, at beam width spacing, after 1 hour of integration. 



be taken to chop in the direction of the faintest extended emission 
(here at a position angle of 20° east-of-north and 2.5' throw); there 
are some minor d iscrepancies with the 350p.m polarization map of 
ISchleuning et al]l2000l) . which can plausibly be explained by their 
larger (and more suitable) 6.3' chop throw. Examination of a JCMT 
archival 850 [im scan-map shows that the faintest regions with plot- 
ted vectors in Figure [5] have in the worst cases up to 25% relative 
flux in one of the off-beams. Although no large-scale polarization 
data are available, we can estimate that a 10% polarization at this 
off-point would seriously contaminate a 2.5% on-source vector. 

Projects on regions such as W3 are producing very valuable 
information on different magnetic morphologies in clouds with 
more or less active star formation and varying core and star masses 
I Matt hews & Wilsonll2002|. | 2000t ICoppin et al]l200fi iDavis et alj 
l200filFiege & Pudritzll200cl) . Also, the marked decrease in polar- 
ization with increasing flux can be used to test models of field lines 
that are tangled or turbulent within the beam, and/or models of 
grains that are less well aligned in dense cores lOstriker et all200lt 
lArce etalll99l . 

The faintest targets detected with jiggle map polarimetry have 
fluxes of about 0.2 Jy/beam and a few percent polarization. This 
covers a wide range of Gala ctic interstellar clou ds and also the 
closest star-forming galaxies iGreavesetalE oOO). Observing ef- 
ficiency compared to previous point-by-point magnetic field map- 
ping has been greatly increased: for example, the 800vim polarime- 
try of the W3-IRS4 core reported by Greaves et ah Nl99St) required 
one hour of observation per spatial point, whereas the data in Fig- 
ure|5|took a similar time for 78 detected points going down to much 
lower flux limits. 



5.3 Polarimetric Scan-mapping 

Scan map polarimetry is still under development, but initial tests 
on the polarized synchrotron emission of the Crab Nebula gave 
very good results (Figure [6j. The 850 \im polarization was very 



similar between scan map observations and a standard jiggle map 
(of the central 3 arcmin); percentages and angles differed by only 
1% and 6° (the latter is probably dominated by rotation of a few 
degrees during each 30 minute map sequenc e). The angles also 
agree very well with 9 mm results published bv lFlett & Henderson! 
1 1979). However, the SCUBA data improve our understanding of 
the (sub)millimetre polarization pattern from 90 to 15 arcsecond 
resolution, greatly clarifying the geometry around features such as 
the emission 'arm' to the north-west. The average flux and polar- 
ization for the Crab Nebula are 1 Jy/beam and 25%, and 5<r 
detections were obtained in one hour of observing. A more typical 
source with an equivalent polarized flux, such as a 5 Jy/beam, 5% 
polarized dust cloud, should be detectable in the same time. 



6 CURRENT AND FUTURE DEVELOPMENTS 

Work is ongoing on automating scan-map polarimetry and expand- 
ing the use of 450 \im polarimetry. The latter will give a factor 
of two improvement in spatial resolution over 850 [im but not all 
the IP data have been obtained. Fast corrections for atmospheric 
chang es are also b eing tested, utilising a water vapour radiometer 
(Wiedner et al. 2001) in the JCMT's receiver cabin. This monitors 
the profile of the 183 GHz atmospheric absorption line about every 
second, and thus gives a very precise transmission measurement 
along the same line-of-sight as the source. Some tests have also 
been made of a 'continuous spinning' mode, as an alternative to 
stepping the waveplate to discrete angles. By analysing the power 
spectrum of the signal, it is possible to detect the intensity of the 
polarized signal (a modulation every 90°), and the position angle 
can also be measured by retaining information on orientation as a 
function of time. This should be a very fast method of detecting 
strong sources, and rapid spinning (e.g. > 1 Hz) would eliminate 
most problems with the sky changing before the complete polarized 
modulation has been detected. 

Finally, a replacement camera is currently being built by the 
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UK Astronomy Technology Centre. With SCUBA-2 and a new po- 
larimeter, there should be a factor of 2.5 decrease in NEFD, a factor 
of 16 increase in imaging speed due to instantaneous full-sampling 
(no jiggling) and a 15-fold increase in field-of-view, combining 
to give an improvement in polarization imaging power for large 
sources of ~ 1500. 



7 SUMMARY 

An efficient polarimeter coupled with a very sensitive submillime- 
tre camera has facilitated a major advance in ground-based imaging 
polarimetry. The variety of science targets has also expanded, with 
almost every category of Galactic source of dust and synchrotron 
emission becoming available, as well as a number of extragalac- 
tic objects. Investigations of the dust polarization spectrum can 
also be made, especially in conjunction with 60-350 \im polar- 
ization data from the Kuiper Airborne Observatory and the CSO 
( Hildebrand et al. 2000; Dotson et al. 2000). This reveals unique 
information about the relation of grain alig nability an d emissiv- 
ity, particularly in mix ed grain populations iGreaves et alJll999l 
Hilde brand et alll999l) . Finally the multi-wavelength capability of 
the SCUBA Polarimeter covers both the high-frequency end of syn- 
chrotron spectra and the long-wavelength tail of dust spectra, where 
both effects are generally optically thin and easiest to interpret. 
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